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Abstract 

Several aspects of the theory of the coexistence of phases and equilibrium 
forms are discussed. In section 1, the problem is studied from the point of 
view of thermodynamics. In section 2, the statistical mechanical theory is 
introduced. We consider, in particular, the description of the microscopic 
interface at low temperatures and the existence of a free energy per unit area 
(surface tension). In section 3, a proof is given of the microscopic validity of 
the Wulff construction in a simplified model of a two-dimensional interface. 
Finally, in section 4, the roughening transition and the formation of facets in 
an equilibrium crystal are studied. Appendices A and B concern, respectively, 
the first and second points of section 2, mentioned above. 
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S. Miracle-Sole / Interfaces 

1 Introduction 

When a fluid is in contact with another fluid, or with a gas, a portion of 
the total free energy of the system is proportional to the area of the surface 
of contact, it is also proportional to a coefficient, the surface tension, which 
is specific for each pair of substances. Equilibrium, which is the condition 
of minimum free energy in the system, will accordingly be obtained by the 
utmost possible reduction of the surfaces in contact. 

We must understand here, not an absolute but a relative minimum, de- 
pending on the material circumstances of the case that is considered. Math- 
ematicians restrict the term minimal surfaces to the former cases, or more 
generally, to all cases in which the mean curvature is zero. Other surfaces 
which are only minimal with respect to the volume contained, are called sur- 
faces of constant mean curvature. From the fact that we may extend a soap 
film across any ring of wire, however the wire be bent, we see that there 
is no end to the variety of the minimal surfaces that can be constructed or 
imagined. To fit a minimal surface to the boundary of any given closed curve 
in the space is a problem formulated by Lagrange, commonly known as the 
problem of Plateau, who solved it with the soap films. 

When we have three fluids in contact with one another, new interesting 
situations arise. Suppose that a drop of some fluid, o, float on another fluid, 
6, while both are exposed to air, c. Then we have three different surfaces 
of contact, and the total free energy of the system consists of three parts, 
associated to these three surfaces. Equilibrium will be reached by contracting 
the surfaces whose surface tension happens to be large and extending those 
where it is small. A drop of fluid a, will exist provided its own two surface 
tensions exceed the surface tension between fluid b and the air, i.e., provided 
that 

Tab + T ac > T bc 

If equality is attained, then a film of fluid a is formed, a situation which is 
known as perfect wetting. For instance, floating on water, turpentine gathers 
into a drop and olive oil spreads out in a film. 

When one of the substances involved is anisotropic, such as a crystal, 
the equilibrium is not determined as it was in the case of a fluid, by the 
condition of minimal area. The contribution to the total free energy of each 
element of area of the crystal surface depends on its orientation with respect 
to the crystalline axes. 

Gibbs' article "On the equilibrium of heterogeneous substances" (Gibbs 
1875) is a theoretical system of thermodynamics, of nearly 300 pages, derived 
from the first and second law. In this work, Gibbs introduced the principle 
of minimum of total surface free energy. He also showed the role of the 
anisotropic surface tension for determination of the shape of a crystal in 
equilibrium. Concerning the formation of facets, he remarks: 
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"the value of r varies greatly with the direction of the surface with 
respect to the axes of crystallization and in such a manner as to 
have sharply defined minima (The differential coefficients of r with 
respect to the direction-cosines of the surface appear to be discon- 
tinuous functions of the latter quantities)." 

He also points out the complexities of actual growth of a crystal and suggests 
that only small crystals have an ideal equilibrium form: 

"it seems not improbable that the form of very minute crystals in 
equilibrium is principally determined by the condition that the free 
energy shall be a minimum for the volume of the crystal, but as they 
grow larger the deposition of new matter on the different surfaces 
will be determined more by the nature (orientation) of the surface 
and less by their size in relation to the surrounding surfaces. As a 
final result, a large crystal thus formed will generally be bounded 
by those surfaces alone on which the deposition of new matter takes 
place least readily" 

Curie (1885) considered also, independently of Gibbs and using a differ- 
ent method, a crystal in equilibrium with its own vapor phase, arriving at 
the same conclusion. For a crystal of N faces, of area Fj and specific surface 
free energies Tj, he found that 

/ = n-Fi + t 2 F 2 + ... + t n F n 

has the minimum value for a given volume. As an application he examined 
two problems, a crystal taking the form of a parallelepiped, and a crystal 
taking the form of a of a cubo-octahedron. The second problem, for which 
Curie gives only the result, is not obvious. 

The article by Curie was the point of departure for Wulff's classical 
experiments on crystal growth. He reported the results in a paper (Wulff 
1901), first published in Russian in 1895. His principal conclusions include 
the celebrated Wulff's theorem on the shape of the equilibrium crystal: 

"The minimum surface energy for a given volume of a polyhedron 
will be achieved if the distances of its faces from one and the same 
point are proportional to their capillary constants." 

The term capillary constants was used for the surface tension. Denoting 
these distances by h i: this means that 

hi _ h 2 _ h N 

Wulff himself supported his principle mainly by its consequences (relating the 
rate of growth in various directions to the corresponding surface tensions). 
His attempt at a general proof failed. He surmised that the surface free 
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energy / (given by Curie's formula) and the volume V are related to hi by 
the following expressions 

f=pJ2nhl V = qJ2hf 

where p and q are constants. However, the simplest example of a paral- 
lelepiped contradicts these statements. 

The complete proof of Wulff theorem was presented by von Laue (1943), 
Herring (1951), and others. 

The surfaces tensions are dependent upon the geometric distribution of 
the particles making up the crystal structure. In this way the Wulff theorem 
establish a relation between the forms and the structure of crystals. Although 
it is not easy to find equilibrium crystals in nature we may think, following 
Schneer (1970), that the ensemble of the forms of a mineral species should 
approach equilibrium. The distribution of these forms in geological time and 
space is the outcome of a world-experiment establishing the probability of 
occurrence of each particular form. 

A complete map of the mineral kingdom known to science was presented 
in Goldschmidt Atlas (1913-1923). Each of the nine volumes of illustrations 
was accompanied by a volume of morphological data. Ironically, the first 
volume of the Atlas appeared in 1913, one year after the discovery of X- 
ray difraction. Its purpose was to provide data for the determination of the 
molecular forces within the crystal from the crystal form. 

It is only in recent times that equilibrium crystals have been produced 
in the laboratory. Most crystals grow under non-equilibrium conditions, as 
predicted by Gibbs, and is a subsequent relaxation of the macroscopic crystal 
that restores the equilibrium. This requires transport of material over long 
distances and the time scales can be very long, even for very small crystals. 
On has been able to study, however, metal crystals in equilibrium of size 
1-10 micron (Heyraud and Metois 1980, 1983). Equilibration times of a few 
days were observed. 

A very interesting phenomenon that can be observed on equilibrium 
crystals is the roughening transition. This transition is characterized by the 
disappearance of a facet of a given orientation from the equilibrium crystal, 
when the temperature attains a certain particular value. Roughening tran- 
sitions have been found experimentally, first, in negative crystals, i.e., vapor 
bubbles included in a crystal, in organic substances (Pavlovsca and Nenov 
1971, 1973). The best observations have been made on helium crystals in 
equilibrium with superfluid helium (Balibar and Castaign 1980, 1985, Keshi- 
shev et al. 1981, Wolf et al. 1983), since the transport of matter and heat 
is extremely fast in the superfluid. Crystals grow to size of 1-5 mm and 
relaxation times vary from milliseconds to minutes. Roughening transitions 
for three different types of facets have been observed. 



4 



S. Miracle-Sole / Interfaces 
1 Thermodynamics of equilibrium forms 

1.1 Variational problem 

Our aim is to find the shape of a droplet of a phase c, the crystal, inside 
a phase m, called the medium, when both phases are in equilibrium. The 
phase c need not be a crystal, this terminology is adopted to stress that 
we consider that c has anisotropic properties. This shape is obtained by 
minimizing the surface free energy between the phases c and m. 

Let n be a unit vector in R d and consider the situation in which the 

phases c and m coexist over a hyperplane perpendicular to n, we denote by 

r(n) the surface tension, or free energy per unit area, of such an interface. We 

consider the surface tension r(n) as a positive function on the unit sphere 
S d-i c R d but 

we do not require now the symmetry r(n) = r(— n). In 
the case under consideration r(n) is obviously a symmetric function, but 
non-symmetric situations are also physically interesting and they appear, for 
instance, in the case of a sessile drop on a wall to be be discussed in section 
1.6. 

We denote by B the set of R d occupied by the phase c, and by dB the 
boundary of B. The total surface free energy of the crystal is given by 

r(dB) = [ r(n(0)ds C (1.1) 

Here n(£) is the exterior unit normal to dB at £ and ds^ is the element of 
area at this point. We suppose that the boundary of B is sufficiently smooth 
so that the normal is defined for almost all s in dB. 

The variational problem is to minimize (1.1) under the constraint that 
the total (Lebesgue) volume \B\ occupied by the phase c is fixed. Given a 
set W, we say that the crystal B has shape W if after a translation and a 
dilation it equals W. The solution of the variational problem, known under 
the name of Wulff construction, is given below. Notice that the problem is 
scale invariant, so that if we can solve it for a given volume of the phase c, 
we get the solution for other volumes by an appropriate scaling. Let W be 
defined by 

W = {x G R d : x- n < r(n)} (1.2) 

where the inequality is assumed for every n e S d_1 (x-n denotes the euclidean 
scalar product). We shall show that this set, which will be called the Wulff 
shape, gives the optimal shape for the crystal. 

Theorem 1.1. Let W be the Wulff shape (1.2) for the surface tension 
function r(n). Let B C R d be any other region with sufficiently smooth 
boundary and the same (Lebesgue) volume as W. Then 

t{8B) > r(dW) (1.3) 
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with equality if and only if B and W have the same shape. 

A proof of this Theorem, based on some version of the isoperimetrical 
inequality, will be given at the end of section 2.2 (in this proof we follow 
Taylor 1987). 

We now state some elementary properties of W. Being defined in (1.2) as 
the intersection of closed half-spaces the Wulff shape W is a closed, bounded 
convex set, i.e., a convex body. 

Among the functions r(n), which through (1.2) define the same shape 
W, there is a unique function having the property that all planes {x e R d : 
x ■ n = r(n)}, associated to all different unit vectors n, are tangent to the 
convex set W. This function is given by 

ryy (n) = sup (x • n) (1.4) 

and is called the support function of the convex body W. If we consider an 
arbitrary function r(n) defining the Wulff shape W, it can be that some of 
these planes do not touch the set W. In other words, we have always 

r(n) > r w (n) (1.5) 

but the inequality sign can hold for some n. The support function ryy(n) is 
the smallest function on the unit sphere which gives the Wulff shape W. 

1.2 Geometric inequalities 

We can state as follows the isoperimetric property of convex bodies. 

Theorem 1.2. Let W C H d be a convex body and r>v(n) the cor- 
responding support function. For any set B C Ti d , with sufficiently smooth 
boundary, we have 

rw(dB) > d \W\ 1,d \B\^ d - 1)/d (1.6) 

where \W\, \B\, denote the (Lebesgue) volumes ofW, B, respectively, and 
Tyv(dB) is the surface functional defined by (1.1). The equality occurs only 
when B and W have the same shape. 

Let us recall the following classical statement, known as the isoperimetric 
property of the circle. 

Proposition. The area F and the length L of any plane domain with 
rectifiable boundary satisfy the inequality 

L 2 > 4nF (1.7) 

The equality sign holds only for the circle. 
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If for the set W we choose D, the closed circle of unit radius with center 
at the origin, then the corresponding support function T£>(n) is equal to the 
constant 1, and we obtain the isoperimetric property. 

Mathematicians were convinced that the isoperimetric property of the 
circle is valid since ancient times. The known proofs of inequality (1.7) may 
be grouped together around several main approaches. Three of them go 
back to Steiner, who gave the first non rigorous proofs of the isoperimetric 
property of the circle on the plane and of the sphere and proved a similar 
property of the ball in space. 

We shall describe here another approach based on the notion of vector 
addition of sets due to Minkowski. This is a very general approach which 
easily leads to the Brunn-Minkowski inequality and, as a particular case, to 
the isoperimetric inequality. 

To every pair of non-empty sets A,5c R d their vector Minkowski sum 
is defined by A + B = {a + b : a G A, b E B}. If A, B are compact sets, then 
A + B is compact. In this case each of the sets necessarily has a volume, its 
Lebesgue measure. The following Theorem is known as the Brunn-Minkowski 
inequality. 

Theorem 1.3. For non empty compact sets A,Bc H d 

lA + Bl 1 ^ > \A\ 1 / d + \B\ x / d (1.8) 

and the equality sign in (1.8) holds only in the following cases: \A + B\ = 0, 
one of the sets A, B consists of a single point, and A and B are two convex 
bodies with the same shape. 

A proof of this Theorem can be found in the book by Burago and Zal- 
galler (1988). First one proves by direct computation that the inequality 
holds in the case in which A and B are parallelepipeds with sides parallel to 
the coordinate axis. The validity of the inequality is then extended by induc- 
tion to all finite unions of such parallelepipeds and, finally, to all compact 
sets by an appropriate limit process. 

The following geometrical remark will next be used to prove the isoperi- 
metric inequality. 

Remark. Let W be a convex body in R d . Given any set B C R d , with 
sufficiently smooth boundary, we can express the functional Tw(dB) as 

\B + \W\-\B\ 
t w (8B) = hm (1.9) 

A->-0 A 

where AW denotes the homothetic set {Ax : x e W}. If B C R d is bounded, 
the functional ryy; is well defined. In particular, (1.9) shows that 

Tw (dW)=d\W\ (1.10) 
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Proof of Theorem 1.2. We first prove equation (1.9). In fact, if H(n) 
denotes the half space whose oriented unit normal is n, we obtain from the 
definition of ryy 

r w (n) = dist (dH(n), d(H(n) + W)) (1.11) 

We may then write 

r w (dB) = ! dist (<9#(n(£)), d(H(n(£)) + W))ds c 
JiedB 

= [ lim(l/A)dist(f,0(£ + \W))dsz (1.12) 

This last expression coincides with the right hand side of (1.9), and proves 
the validity of this equation. 

Then inequality (1.6) follows by applying the Brunn-Minkowski inequal- 
ity to \B + XW\. This gives 

\B + XW\ - \B\ > (\B\ 1 ' d + X\W\ 1/d ) d -\B\>d\ \B\^ d -^/ d \W\ 1/d (1.13) 

which, taking equation (1.9) into account, ends the proof of inequality (1.6). 

By taking limits we can deduce (1.6) from (1.8). However, conditions 
for the equality to occur in (1.6) cannot be found in this manner. For a proof 
that W is the unique minimum see Burago and Zalgaller (1988). 

Proof of Theorem 1.1. The proof can be obtained as a consequence of 
the Theorem 1.2. Let W be the Wulff shape corresponding to the function 
r. Then 

r(dB) > t w (8B) > d |W| 1/d |B| (d - 1)/d (1.14) 

taking into account equation (1.9) and the isoperimetric inequality (1.6). 
But, when B = W, we have 

T (&W) = r w (dW) = d \W\ (1.15) 

The first equality in (1.15) follows from the fact that rw(n) ^ T(n) only for 
the unit vectors n for which the planes x • n = r(n) are not tangent to the 
convex set W. The second equality follows from equation (1.9). Therefore 

r(dB) > r(dW) (\B\/\W\) {d - 1)/d (1.16) 

which, when \B\ = \W\, gives the stated inequality (1.3). The equality in 
(1.3) corresponds to the equality in (1.6). The Theorem is proved. 
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1.3 Convexity properties 

It is expected, from thermodynamic reasons, that the surface tension 
r(n) is equal to the support function of the associated Wulff shape W. We 
next comment on this point, first discussed by Dobrushin and Shlosman 
(1992) and Messager et al. (1992). 

Let A , . . . , Ad e H d be any set of d + 1 points in general position and, 
for i = 0, . . . , d, let Aj be the (d — 1) -dimensional simplex defined by all these 
points except Ai. Denote by |Aj| the (d— 1) -dimensional area of Aj and by 
rij the unit vector orthogonal to Aj. The first vector, no, is oriented toward 
the exterior of the simplex Ao, . . . , Ad, while the others, rij, i = 1, . . . , d, are 
oriented inside. We say that r(n) satisfies the pyramidal inequality if 

d 

|A |r(n )<^|A l |r(n i ) (1.17) 

i=i 

for any set Ao, . . . , Ad- We introduce the function on R d defined by 

/(x) = |x|r(x/|x|) (1.18) 

If this function is convex, it satisfies 

/(ax) = a/(x) (1.19) 
/(x + y) </(x) + /(y) (1.20) 

for any a > and any x and y in R d . The first condition, that /(x) 
is a positively homogeneous convex function of degree 1, is clear from the 
definition. The second is equivalent in this case to the convexity condition. 

Theorem 1.4. The following propositions are equivalent 

1) r(n) satisfies the pyramidal inequality 

2) the function /(x) is convex 

3) r(n) is the support function of the convex body W 

Proof. For simplicity we shall restrict ourselves to the three dimensional 
case. Assume that condition 3) is satisfied and let us consider the set W and 
its tangent planes orthogonal to no, . . . , 113. These planes form a pyramid 
whose four vertices are denoted by Aq, A±, A2, A3. The notation is chosen in 
such a way that the plane A±,A2, A3 separates the origin O from the vertex 
Aq. This is possible, for all such pyramids, if and only if W is convex. We 
consider also the four pyramids having as vertices the origin and three of 
these points. The volumes of these pyramids satisfy 

|OAiA 2 A 3 | < \OA A 2 A 3 \ + \OA A 3 A 1 \ + |OA AiA 2 | (1.21) 

Taking into account that the heights of these pyramids are r(n ), T(ni), 
r(n 2 ) and r(n 3 ), respectively, we get the pyramidal inequality for the pyra- 
mid Ao, Ai, A2, A3. This shows that conditions 1) and 3) are equivalent. 
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Next, we observe that 

|A | n = |Ai| ni + |A 2 | n 2 + |A 3 | n 3 (1.22) 

Thus the pyramidal inequality says that 

/( Xl + x 2 + x 3 ) < /( Xl ) + /(x 2 ) + /(x 3 ) (1.23) 

with x/j = |Afe| rifc, k = 1, . . . , 3. Then, to prove (1.20), it is enough to find 
pyramids such that 

|Ai|ni^x, |A 2 |n 2 ^y, |A 3 |n 3 ^0 (1.24) 

for all x and y. This can easily be done and proves the equivalence of 1) and 
2). 

The pyramidal inequality may be interpreted as a thermodynamic sta- 
bility condition and thus also the convexity of /(x). If one supposes that 
|Ao| t(iio) is greater than the right hand side of inequality (1.17) this would 
make the interface bounded by the sides of A unstable and difficult to real- 
ize. 

Let us mention that the term of support function is usually attributed to 
the convex function /(x). In fact a support function can be defined, by the 
same expression, for any convex set in R d , and it is a real (or +oo) valued 
function, positively homogeneous and convex. For a non-empty convex body 
the convex function /(x) is finite and hence (see e.g. Ruelle 1969, Proposition 
3.3.4) Lipschitz continuous. 

If /(x) = /(— x) the convex body W is symmetric with respect to the 
origin. Now, if W is a symmetric closed bounded convex set and the origin 
belongs to the interior of W, then the support function of W is finite every- 
where, symmetric and strictly positive except at the origin, i.e., it is a norm. 
Let then /(x) = sup y (x ■ y//(y)) be the dual norm. In this case, the Wulff 
equilibrium shape W can be interpreted as the unit ball with respect to the 
dual norm 

W = {x G R d : /(x) < 1} (1.25) 

1.4 Legendre transformation 

The Wulff construction can also be viewed as a geometrical version of 
the Legendre transformation. We consider the function /(x), defined from 
the surface tension r(n) by equation (1.18), but we do not require now that 
/ is the support function of a convex body. From definition (1.2), we get 

W = {xGR d :/*(x) <0} (1.26) 
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where /* is the Legendre transform of /, 

/*(x)=sup(x.y-/(y)) (1.27) 
y 

Actually /*(x) = if x e W, and /*(x) = oo otherwise. The support 
function of W is given by /**, the Legendre transform of /*. 

This point of view was developed by Andreev (1981). For simplicity we 
consider the case d = 3 and assume also that W is symmetric with respect 
to the origin. We introduce a function <p on R 2 such that the graph of 
xs = <p(x\,x 2 ), for xs > 0, coincides with the boundary <9W of the crystal 
shape. Since W is a convex body, <p is a concave function, and 

W = {x G R 3 : -(p(-x u -x 2 ) <x 3 < <p(x 1 , x 2 )} (1.28) 

In the present context this means that ~f is the Legendre transform of the 
projected surface tension t p = (l/n^r, considered as a function on R 2 of 
the slopes x\ = ni/ri3,x 2 = n 2 /n 3 ). In other words, 

r p (x 1 ,x 2 ) = f(xi,x 2 ,l) I 1 - 29 ) 
Indeed, from equations (1.26) and (1.27), we see that 

-(f(xi, x 2 ) = sup (xiyi + x 2 y 2 - r p (yi, y 2 )) (1.30) 

Formula (1.28) is known as the Andreev construction. The interest of 
this approach comes from the fact that cp, and hence, the crystal shape 
itself, may be regarded as the free energy associated to a certain statistical 
mechanical Gibbs ensemble. We shall consider this ensemble in section 3. 

1.5 Facets in the equilibrium crystal 

Another consequence of the convexity properties, which will next be 
discussed (following Miracle-Sole 1995a), concerns the formation of facets 
in the equilibrium crystal. The facets of a crystal have certain particular 
orientations. Let no be the normal corresponding to one of these orientations, 
and denote by t(6, (p) the function r(n) expressed in terms of the spherical 
co-ordinates, < 6 < < (p < 2n, of n, the vector n being taken as the 
polar axis. 

Theorem 1.5. Assume that the convexity condition is satisfied. A 
facet orthogonal to the direction no appears in the Wulff shape if, and only 
if, the derivative dr(6, <p)/d9 is discontinuous at the point 9 = 0, for all (j>. 
Moreover, the one-sided derivatives Or (6, <p)/d9 , at 9 = + and 9 = 0~ , 
exist, and determine the shape of the facet. 
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Proof. In terms of the function / defined by (1.18), the Wulff shape W 
is the set of all x = (xi, x 2 , x 3 ) G R 3 such that 

xiyi + x 2 y2 + x 3 y 3 < f(y) (1.31) 

for every y G R 3 . If the coordinate axes are placed in such a way that 
no = (0,0, 1), the plane x 3 = r(0) (where r(0) is the value of r for 9 = 0) 
is a tangent plane to W. The facet T is the portion of this plane contained 
in W. These facts follow from the convexity of /, which implies that / is 
the support function of W. According to (1.31), the facet J 7 consists of the 
points (xi, X2, t(0)) G R 3 such that 

xiyi + x 2 y2 < f(yi, 2/2, 2/3) - 2/3t(0) = f(y 1: y 2 , 2/3) - 2/3/(0, 0, 1) (1.32) 

for all y = (2/1,2/2, 2/3)- Or, equivalently, such that 

xiyi + x 2 y2 < £(2/1,2/2) = inf (/(2/1, 2/2, 2/3) -2/3/(0,0, 1)) (1.33) 

2/3 

Restricting the infimum to 2/3 = 1/A > 0, and using the positive homogeneity 
and the convexity of /, one obtains 

g(y u 2/2) = ^ lim (1/A) (f(X yi , Xy 2 , 1) - /(0, 0, 1)) (1.34) 

This implies that g is a positively homogeneous convex function on R 2 . De- 
fine 

fj,{4>) = g(cos(f), sin0) (1.35) 
From (1.34), and taking A = tan 9, one gets 

fj,(4>) = lim (l/sin6>) (/(sin^cos^, sin6>sin^, cos^) — cos6'/(0, 0, 1)) 
= lim (1/9) [r(9, 0) - r (0)) = (d/d6) e=a +T(0, 4>) (1.36) 

Similarly 

M0 + tt) = £(-cos0,-sin0) = -(d/de) e=0 -T(e t (f>) (1.37) 

Equations (1.33) and (1.35) show that both one-sided derivatives of exit and, 
from the convexity of g, it follows that 

M0 + 7r)<M0) (1.38) 

Thus, the hypothesis of the discontinuity of the derivative dr/d9, at 9 = 0, 
implies the strict inequality in (1.38) and shows that the convex set J 7 , given 
by (1.33), has a non-empty interior. This ends the proof of the Theorem. 

12 



S. Miracle-Sole / Interfaces 
Notice that, according to condition (1.33), the shape of the facet is given 

by 

T = {x G R 2 : xm < //(m)} (1.39) 

where the inequality is assumed for every unit vector m = (cos0, sin0). 

If r(9,(j)) = r(n — 6,<p) the crystal shape is reflection symmetric with 
respect to the plane xs = 0. Then //(m) = m) (assuming also the 
symmetry with respect to the origin), and the facet J 7 has a center at the 
point Po = (0,0, r(0)). This is not always the case, however, and the point 
Po can also be outside J 7 . 

1.6 Crystal on a wall 

The variational problem, considered in section 1.1, can be solved by 
purely geometrical means also in the presence of walls. We refer the reader 
the to paper by Kotecky and Pfister (1994) for a general analysis of these 
situations, and consider here, as a basic illustration, the case of a crystal on 
a single wall (Winterbottom 1967). 

In this case, the surface tension of the boundary of the phase c in contact 
with the medium differs from that arising with the wall w. We suppose that 
the wall is described by a plane perpendicular to the unit vector n , and that 

E = {x G R d : x-n < a} (1.40) 

is the half-space where we have the phases c and m. The surface tension 
between the phases c and m is denoted by r(n). The relevant physical 
quantity in this problem is the difference 

a(n ) = t cw (yi ) - T mw (n ) (1.41) 

where t cw and r mw are the surface free energies of the phase c with the wall 
and of the phase c with the wall, respectively. Since cx(no) is a difference of 
free energies, it may be positive or negative. 

When a(n ) > 7"(n ) we have a drying situation: in equilibrium, it is 
preferable that the medium m occupies the place between the wall and the 
phase c, and consequently the phase c is not in contact with the wall. On the 
other hand, when cr(no) < —t(yiq) we have a (complete) wetting situation: 
in equilibrium, the phase c forms a layer between the wall and the medium 
m. 

In all other cases we speak of partial drying or partial wetting. Then, if 
B is the shape of the phase c, the total surface free energy is given by 

t{8B) = \dB n dE\ <T(n ) + / r(n(f))ds € (1-42) 

where the first term is the contribution of the boundary of the crystal along 
the wall and the second is the integral over the remaining part of the bound- 
ary of the crystal. The solution of the variational problem for the functional 
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(1.42) is the Winterbottom shape W, and can be obtained as follows. One 
first constructs the Wulff shape W r associated to the function r(n) which 
corresponds to the ideal shape of the free crystal (formula (1.2)). Then we 
take the intersection of this set with the half space defined by the plane 
orthogonal to no, at a distance cr(no) from the origin, i.e., 

W = W T n{xGR d :x-no< ( j(no)} (1.43) 

We notice that the set W is a convex body, still defined according to formula 
(1.2), but with the surface tension being replaced by a(n ) when n is the 
oriented normal to the interface. Thus, equation (1.43) follows, by arguing 
as in the proof of Theorem 1.1, from the isoperimetric inequality. 

1.7 Stability of the Wulff solution 

Let us finally comment on the stability of the Wulff shape as the solution 
to the variational problem of section 1.1. In the two-dimensional case, in 
addition to the isoperimetric inequality, the stability of the minimum in 
Theorem 1.1 can be controlled. The following classical inequality, known as 
the Bonnensen inequality, gives an estimate from below of the isoperimetric 
defect in the case of the circle. Suppose that a closed curve of length L 
bounds a plane domain £?, of area F, and suppose that r and R are the radii 
of the incircle and the circumcircle, i.e., the largest circle contained in B and 
the smallest circle containing B. Then 

L 2 -4nF > n 2 (R 2 - r 2 ) (1.44) 

In the general case, when W is a convex body, a similar inequality for the 
isoperimetric defect 

t w (8B) 2 -4 \W\ \B\ (1.45) 

has been proved by Dobrushin et al. (1992). 

There are no estimates of this kind in the three dimensional case. In two 
dimensions we have that if Tw(dB) — Ty^(dW) tends to zero, and \B\ = |W|, 
then the Hausdorff distance 

dist H (dW,B') = max{ sup dist(x,dB), sup dist(x,dW)} (1.46) 

xEdW xedB 

between the boundaries of W and of a suitably chosen translated B' of B, 
tends also to zero. This is true in three dimensions if both W and B are 
convex bodies, because of Diskant inequalities (see for instance Burago and 
Zalgaller 1988), and some appropriate distance could also be suggested be- 
tween general W and B, in order to have the stability of the Wulff solution. 
However, the smallness of the difference Tw(dB) — r>v(9VV) does not im- 
ply anymore the smallness of the Hausdorff distance between dB' and <9W. 
This is because a "hair" attached to the Wulff surface does not contribute 
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significantly to the total surface free energy. From the point of view of ther- 
modynamics, the Wulff shapes with long hairs appear to have almost the 
same free energy as the equilibrium shape. 



2 Interfaces in statistical mechanics 

In this and the following sections we study some aspects of the theory 
of the coexistence of phases in statistical mechanics. 

In the first approximation one can model the interatomic forces in a 
crystal by a lattice gas. These systems consist of lattice cells which may 
be either empty or occupied by a single particle. In a typical two-phase 
equilibrium state there is a dense component, which can be identified as 
the crystal phase, and a dilute phase, which can be identified as the vapor 
phase. The underlying lattice structure implies that the crystal phase is 
anisotropic, while this assumption, though unrealistic for the vapor phase, 
should be immaterial for the description of the crystal-vapor interface. As 
an illustrative example of such systems, the ferromagnetic Ising model will 
be considered. 

The Ising model is defined on the <i-dimensional cubic lattice C = Z d , 
with configuration space O = { — 1, The value a(i) is the spin at the site 
i. The occupation numbers n{i) = (l/2)(<r(i) + 1, which take the values or 
1, give the lattice gas version of this model. The energy of a configuration 
C"A = £ A}, in a finite box A C £, under the boundary conditions 

a G O, is 

H A (a A \a) = - ff (0*0') 

<i,i)nA#0 

where are pairs of nearest neighbour sites and a(i) = a(i) if i A. The 
partition function, at the inverse temperature (5 = 1/kT, is given by 

Z*(A) = ^exp ( - PH A {a A I *)) (2.2) 
The following limit, which is independent of the boundary conditions, 

^-^-m lnZ ' w (2 ' 3) 

defines the free energy per unit volume. 

It is known that this model presents, at low temperatures T < T c , where 
T c is the critical temperature, two distinct thermodynamic pure phases, a 
positively and a negatively magnetized phase (or a dense and a dilute phase 
in the lattice gas language). This means two extremal translation invariant 
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Gibbs states, which correspond to the limits, when A — > oo, of the finite 
volume Gibbs measures 

^(A)- 1 exp(-if A ((T A | a)) (2.4) 

with boundary conditions a respectively equal to the ground configurations 
(+) and (— ), such that a(i) = 1 and a(i) = — 1, for all i G C. On the other 
side, if T > T c , then the Gibbs state is unique (see, for instance, Ruelle 1969, 
Gallavotti 1972b, Miracle-Sole 1976). 

Each configuration inside A can be geometrically described by specifying 
the set Peierls contours, which indicate the boundaries between the regions 
of spin 1 and the regions of spin —1. Unit square surfaces are placed midway 
between nearest-neighbour pairs of sites i and j, and perpendicular to these 
bonds, if a(i)a(j) = —1. The connected components of this set are the 
Peierls contours. Under the above boundary conditions (+) and (— ), they 
form a set of closed polyhedrical surfaces. The contours can be viewed as 
defects, or excitations, with respect to the ground states of the system (the 
constant configurations 1 and —1), and are a basic tool for the investigation 
of the model at low temperatures. 

2.1 The surface tension 

In order to study the interface between the two pure phases one needs 
to construct a state describing the coexistence of these phases. Let A be 
a parallelepiped of sides Li, . . . , Ld, parallel to the axes, and centered at 
the origin of £, and let n = (rii, . . . , nd) be a unit vector in R d , such that 
nd 7^ 0. Introduce the mixed boundary conditions (±, n), for which a(i) = 1 
if i ■ n > 0, and a(i) = — 1 if i ■ n < 0. These boundary conditions force 
the system to produce a defect going transversely through the box A, a big 
Peierls contour that can be interpreted as a microscopic interface. The other 
defects that appear above and below the interface can be described by closed 
contours inside the pure phases. 

The free energy, per unit area, due to the presence of the interface, is 
the surface tension. It can be defined by 

r(n)= lim lim - d In \> (2.5) 

Li,...,L d _!->-oo L d ^-oo pL\...Ld-\ Z^){A) 

Notice that in this expression the volume contributions proportional to the 
free energy of the coexisting phases, as well as the boundary effects, cancel, 
and only the contributions to the free energy of the interface are left. 

Theorem 2.1 The thermodynamic limit r(n), of the interfacial free 
energy per unit area, exists, and is a non negative bounded function ofn. Its 
extension by positive homogeneity, /(x) = |x| r(x/|x|), is a convex function 
on K d . 
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A proof of these statements has been given by Messager et al. (1992) 
using correlation inequalities (this being the reason for their validity for all 
(3). In fact the validity of this Theorem has been proved for a large class of 
lattice systems. 

Moreover, for the Ising model we know, from Bricmont et al. (1980), 
Lebowitz and Pfister (1981) and the convexity condition, that r(n) is strictly 
positive for T < T c and that it vanishes if T > T c . 

2.2 Interfaces at low temperatures 

In Appendix A we give a resume of the theory of cluster expansions. We 
are going to apply this theory to study the low temperature properties of the 
system. Let us first describe the cluster theory of a pure phase. We have seen 
above that a configuration in the box A, with (+) boundary conditions, is 
characterized as a set X = {71, . . . , 7 n } of closed mutually disjoint contours. 
It is easy to see that, up to a constant, the energy of the configuration is 
equal to twice the total area of the contours. Besides the mutual exclusion 
there is no other interaction between the contours. If to each contour we 
attribute the weight 0(7) = exp(— 2/3|7|), where I7I is the area of the contour 
7, we get, for the partition function Z^ + \h), an expression like equation 
(A.l), in Appendix A. The system can be viewed in terms of contours as a 
system of polymers, with weight ^(7), and the convention that two contours 
are compatible if they do not intersect. 

Let us next consider the convergence condition, Theorem A. 2 in Ap- 
pendix A. We introduce ^(7) = exp(— &I7I) and remark that the number of 
contours cutting a given lattice bond and of area I7I = n is bounded by 
K n , where K is certain constant (one can take K = 3 using appropriate 
definitions). Then, we have 

7' 7^7 n 

This shows that the convergence condition (A. 6) is satisfied if the temper- 
ature is sufficiently low, i.e., there exist a constant /3 such that (A. 6) is 
satisfied if /3 > (3 . Using formula (A. 4) and estimate (A. 7), one can then 
derive a convergent series expansion for the free energy of the system. Also, 
from (A.4) and (A. 7), it follows that 

Z (±,n) (A)/z (+) (A) = ^ e -2/3|i| exp (_ J- a T (X)0 x ) (2.7) 

where the first sum runs over all microscopic interfaces X, compatible with 
the boundary conditions. The sum in the exponential runs over all clusters 
of contours that intersect X. Each term in the first sum gives a weight 
proportional to the probability of the microscopic interface. From this, a 
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representation of the microscopic interface itself as a polymer system, can be 
derived. 

Consider now the three-dimensional Ising model with the (±, no) bound- 
ary conditions, associated to the direction no = (0,0, 1), which define inside 
the box A a horizontal interface. At low temperatures T > 0, we expect 
the microscopic interface corresponding to these boundary conditions, which 
at T = coincides with the plane is = —1/2, to be modified by small de- 
formations. The microscopic interface can then be described by means of 
its defects, or excitations, with respect to the interface at T = 0. These 
defects, called walls, form the boundaries (which may have some width), be- 
tween the smooth plane portions of the interface. In this way the interface 
structure, with its probability distribution defined by equation (2.7), may 
be interpreted as a "gas of walls" on a two-dimensional lattice. There is 
an interaction between the walls, coming from the fact that the interface is 
surrounded by the (+) phase from above and the (— ) phase from below. But 
the exponential function in equation (2.7), gives a mathematical description 
of this interaction in terms of clusters of contours. This allows us to define 
a "gas of aggregates of walls and clusters" as a polymers system on the two- 
dimensional lattice (see the original articles quoted below and, for instance, 
Miracle- Sole 1995a, Appendix, for a short review). 

Dobrushin (1972) proved the dilute character of this gas at low temper- 
atures, which means that the interface is essentially flat (or rigid). The con- 
sidered boundary conditions yield indeed a non translation invariant Gibbs 
state. Furthermore, cluster expansion techniques have been applied by Bric- 
mont et al. (1979), to study the interface structure in this case (see also 
Holicky et al. 1988). 

The same analysis applied to the two-dimensional model shows a dif- 
ferent behavior at low temperatures. In this case the walls belong to a 
one-dimensional lattice, and Gallavotti (1972) proved that the microscopic 
interface undergoes large fluctuations of order \fL[. The interface does not 
survive in the thermodynamic limit, A — >■ oo, and the corresponding Gibbs 
state is translation invariant. Moreover, the interface structure can be stud- 
ied by means of a cluster expansion for any orientation of the interface (see 
also Bricmont et al. 1981). 

In the three-dimensional case, the description of the microscopic inter- 
face for any orientation n, leads to very difficult problems of random surfaces. 
It is possible, however, by introducing a new polymer system, to analyze also 
the interfaces which are near to the particular orientations n , discussed 
above. This problem, related to the formation of facets in the equilibrium 
crystal, will be treated in section 4. 
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3 On the microscopic proof of the Wulff construction 

In this section we consider a simplified model of the two dimensional 
interface. We represent the interface as a single-valued function over a ref- 
erence line. At each site i of the lattice Z an integer variable hi is as- 
signed which indicates the height of the interface at this site. One obtains 
a statistical mechanical model by assigning an energy to each configuration 
h = {ho, hi, /ijv}, in the box < i < N, of length N. The standard 
example is 

N 

H N (h) = J2( 1 + \hi-h i -i\) (3.1) 
i=i 

which corresponds to the solid-on-solid (SOS) model. Then the energy of a 
configuration coincides with the length of the interface. The weight of a given 
configuration, at the inverse temperature (3, is proportional to the Boltzmann 
factor exp ( — /3ii?A (</>)) • Provided that the energy remains unchanged under 
the global shift hi to hi + b, for all i and all b G Z, there is considerable 
freedom of how to choose the energy function. The results below apply, also, 
to more general Hamiltonians of this form. 

The model provides an approximate description of an interface separat- 
ing two phases at equilibrium, such as the positively and negatively magne- 
tized phases of the three-dimensional Ising model. Actually the SOS model 
may be obtained as the limit of the anisotropic Ising model, when we let the 
coupling constant, in the vertical direction, tend to infinity. 

We are going to study the statistical mechanics of the SOS model, with 
the global constraint of having a specified area between the interface and the 
horizontal axis, and to show the existence of the thermodynamic limits and 
the equivalence of the Gibbs ensembles associated with this problem. As a 
corollary, it can be seen that the configurations of a large system follow a well 
defined mean profile with probability one, the shape of the phase boundary at 
equilibrium. This gives a simple alternative microscopic proof of the validity 
of the Wulff construction for such models, first established by DeConink, et 
al. (1989). The fluctuations of the interface around the mean profile have 
recently been studied by Dobrushin and Hryniv (1996) following the same 
approach. 

3.1 Surface tension and Wulff shape 

We introduce first the Gibbs ensemble which consists of all configura- 
tions, in the box of length N, with specified boundary conditions ho = and 
/itv = Y . The associated partition function is given by 

Z 1 (N,Y) = J2e-^ H{h) S(h )S(h N -Y) (3.2) 

h 
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where the sum runs over all configurations in the box and S(t) is the discrete 
Dirac delta (5(t) = 1 if t = and S(t) = otherwise). We define the 
corresponding free energy per site as the limit 

r p (y) = lim --L\nZ 1 {N,yN) (3.3) 

where y = — tan6>, the slope of the interface, is a real number. This free 
energy is called the projected surface tension. The surface tension, which 
represents the interfacial free energy per unit length of the mean interface, is 

r(0) = cos# r p (-taii0) (3.4) 

We introduce a second Gibbs ensemble, conjugate to the previous en- 
semble, whose partition function, in the box of length N, is given by 

Z 2 (N, x) = J2 e- pH{li) e f3xhN S(h ) (3.5) 

h 

where x G R, replaces as a thermodynamic parameter the slope y. We define 
the associated free energy as 

rt x) = ^-jN lnZ2{N ' x) (3 - 6) 



Theorem 3.1. Limits (3.3) and (3.6), which define the above free en- 
ergies, exist. The first, r p , is a convex even function ofy. The second, <p, 
is a concave even function of x. Moreover, r p and —<p are conjugate convex 
functions, i.e., they are related by the Legendre transformations 

-(p(x) = sup (xy - T p (y)) 

y 

T p (y) = sup (xy + <p(x)) (3.7) 

X 



The validity of the above statements follows from the subadditivlty prop- 
erty 

lnZ 1 (N 1 +N 2 ,y(N 1 + N 2 )) > In Zi (AT X , yNi) + In Z\ (N 2 , yN 2 ) (3.8) 

which is straightforward, and standard arguments in the theory of the ther- 
modynamic limit (Ruelle 1969, Galgani et al. 1971). See also, for instance, 
Messager et al. (1992) for a proof of these statements in a more general 
setting. 

The convexity of t p is equivalent to the fact that the surface tension r 
satisfies the pyramidal inequality (see section 1.3). Relations (3.7) between 
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the free energies express the thermodynamic equivalence of the two ensem- 
bles (3.2) and (3.5). These relations imply that the curve z = ip(x) gives, 
according to the Wulff construction, or the equivalent Andreev construction, 
the equilibrium shape of the crystal associated to our system (see section 
1.4). 

The function <p(x) defined by (3.6) is easily computed by summing a 
geometrical series. One introduces the difference variables 

Hi = - hi (3.9) 

for i = 1, N, so that the partition function factorizes and one obtains 

ip(x) = l-/3- 1 ln^e- /3|n|+/3xn (3.10) 

n 

The explicit form of this function is 

if — 1 < x < 1, and <p(x) = — oo otherwise. 

3.2 The volume constraint 

We next define two new Gibbs ensembles. In the first of these ensembles 
we consider the configurations such that h^ = 0, which have a specified 
height at the origin ho = M and which have a specified volume (area) V 
between the interface and the horizontal axis, this volume being counted 
negatively for negative heights, 

N 

V = V(h) = J2K (3.12) 

i=0 

The corresponding partition function is 

Z 3 (N, V,M) = J2 e _/3H(h) 5{h N )8{V{h) - V)5(h - M) (3.13) 

h 

where V and M are understood as their integer parts when they do not 
belong to Z. The second ensemble is the conjugate ensemble of (3.13). Its 
partition function is given by 

Z 4 (N, u^) = J2 e~ PH{h) e^ v W^ h ° S(h N ) (3.14) 

h 

where u G R and |j6R are the conjugate variables. Our next step will be 
to prove the existence of the thermodynamic limit for these ensembles and 
their equivalence in this limit. 
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Theorem 3.2. The following limits exist 

fo{v,m)= lim In Z 3 (N,vN 2 ,mN) (3.15) 

N— too piV 

fa(u,fi) = lim \nZ 4 (N,u,fj,) (3.16) 

and define the free energies per site associated to the considered ensembles. 
Moreover, tfj 3 and —tfj4 are conjugate convex functions 

— ip4,(u, //) = sup (uv + (Um — ip 3 (v, to)) 

t>,m 

^(f , to) = sup (m> + /WTO + ^(-u, /i)) (3-17) 

u,n 

The crucial observation for proving this Theorem is the subadditivity 
property given below. Then we adapt known arguments (Ruelle 1969, Gal- 
gani et al. 1971) in the theory of the thermodynamic limit. For a more 
detailed proof see Miracle-Sole and Ruiz (1994) or Dobrushin and Hryniv 
(1996). 

Proposition. The partition function Z 3 satisfies the subadditivity 
property 

Z 3 (2N, 2(V' + V"),M' + M") > 

Z 3 (N, V, M') Z 3 {N, V", M") e -2/3|M"|/(2iv-i) (3 _ 18) 

Proof. In order to prove this property we associate a configuration h 
of the first system in the box of length 2N ', to a pair of configurations h' and 
h" of the system in a box of length N, as follows 

h2i = K + h'[ , i = 0, . . . , N 
h 2i _ 1 = K_ 1 + h'( ,i = l,...,N (3.19) 

Then h 2N = h' N + h% = 0, h = h' + h' Q ' = M' + M" and 

N N N 

i=l i=0 i=l 

= 2 (V(W) + V(h")) - M" (3.20) 

This shows that the configuration h belongs to the partition function in the 
right hand side of (3.18). Since 

H 2N (h) = H N (W) + H N (h") (3.21) 
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because n<n = n- and nn-\ = n'(, as follows from (3.19), we get 

Z 3 {N,V',M') Z 3 (N,V",M") < Z 3 (2N,2(V + V") - M" , W + M") (3.22) 

Then we use the change of variables 

hi = hi + (M"/(2N - 1)), i = l,..., 2N — 1, 

ho = h , h 2 N = h 2 N = (3.23) 

which gives 

Z 3 (2N, V - M", M) < e 2/3|M"|/(2Ar-i) Z ^2N, V, M) (3.24) 
to conclude the proof. 

3.3 The Wulff construction 

The free energies ip 3 and ip4 , associated to the system with the volume 
constraint, can be expressed in terms of the functions cp and r p introduced 
in section 3.2. We consider first the function ip^. In terms of the difference 
variables (3.9), we have 

N N 

v(h) = Y,hi = J2 m * ( 3 - 25 ) 

i=0 i=l 

and, therefore, 

N 

Z 4 (N,U,ll) = Y[ ^ e -/3|^+/3(«/JV)m i+ /^n^ (g^g) 
i=l rii 

Taking expression (3.10) into account it follows 

JV 

Z 4 (N, u, n) = exp ( - p + /*)) (3-27) 

i=i 

and 

1 N 1 N 

Mu, /*) = ^lim - <p(%i + H)= ^Ijm - £ + **) (3-28) 

i=l i=l 

which, writing the Riemann sum as an integral, implies 

1 f u 

ipAu, fi) = — / (fix + mdx (3.29) 
u J 
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Next we consider the free energy ^3. This function is determined by 
the Legendre transform (3.17). The supremum over u, \i is obtained for 
the values uo,Ho, for which the partial derivatives of the right hand side of 
equation (3.17) are zero: 

v + {dip A /du){u , no) = 

m + {d^ A /dn) {u , no) = (3.30) 

That is, for uo,ho which satisfy 

1 f u ° 1 

—x / <p(x + Ho)dx tp(Ho + uo)=v (3.31) 

Uq Jo u 

— ((p(po) - y>(po + u )) = m (3.32) 

Uq 

Then, from (3.17), (3.28), (3.31) and (3.32), we get 

if> 3 (v, m) = 2^ 4 (w , A*o) - — f(A»o + uo)v{ho + i*o) - Ho<fM) (3-33) 

Uq 

But, using relation (3.7) under the form 

(p(x) = xip'(x) + T p (cp'(x)) (3.34) 
in (3.29) and integrating by parts, we get 

2-04(^0, A*o) = — / T p (<p'(x))dx H ((fio + u )(p(Ho + u ) - ho<p(^o)) 

(3.35) 

Finally, this equation together with (3.33), implies 



ip s (v,m) = — / T p (ip'(x))dx (3.36) 



To interpret these relations, let us consider the graph of the function 

z = <p(x) and an arc BC of this curve between the points B and C having 

the abcises Ho an d Ho + u 0i respectively. We draw the vertical line x = Ho 

passing through the point B, and the horizontal line z = 4>(ho + u o) passing 

through the point C. Let A be the point where these two lines intersect. We 

observe that the right hand side of (3.31) represents the area ABC divided 
2 

by AC . Therefore, the values uq, Ho, which solve (3.31) and (3.32), are 
obtained when this area is equal to v, with the condition coming from (3.32), 
that the slope ABj AC, is equal to m. Then, according to (3.36), the free 
energy i/js(v, m) is equal to the integral of the surface tension along the arc 
BC, of the curve z = <p(x), divided by the same scaling factor AC = uq. 
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We conclude that, for large N, the configurations of the SOS model, 
with a prescribed area vN 2 , follow a well defined mean profile, that is, the 
macroscopic profile given by the Wulff construction, with very small fluctu- 
ations. This follows from the fact that the probability of the configurations 
which deviate macroscopically from the mean profile is zero in the thermo- 
dynamic limit. The free energy associated to the configurations which satisfy 
the conditions above, and moreover, are constrained to pass through a given 
point not belonging to the mean profile, can be computed with the help of 
equation (3.13). Indeed, if the coordinates of the point are (xi,mi), then 
the partition function associated to such set of configurations is given by the 
product 

Z 3 (x 1 N,V',(m-m l )N) Z 3 ((l - x x )N, V", (m - m^N) (3.37) 

with moreover the condition V + V" = vN 2 . As a consequence of equations 
(3.15) and (3.36), and knowing from the macroscopic variational problem 
that the arc BC is the shape which minimizes the integral of t p , we see that 
the probability of this set of configurations tends exponentially to zero as 
N ->■ oo. 

3.4 Remark 

The mathematically rigorous justification of the Wulff construction, in 
the case of the two-dimensional Ising model at low temperatures, is due to 
Dobrushin et al. (1992, see also Dobrushin et al. 1993, for the statement of the 
main result and some ideas of the proof, and Pfister 1991, for another version 
of the same proof). Their results show that, in the canonical ensemble, where 
the total number of particles (or the total magnetization in the language of 
spin systems) is fixed, a (unique) droplet of the dense phase, immersed in the 
dilute phase, is formed. Its shape, when properly rescaled, obeys the Wulff 
principle. Ioffe (1994, 1995) extended these results, in the case of the two- 
dimensional Ising model, to all temperatures below the critical temperature. 

4 Roughening transition and equilibrium shapes 

Let us come back to the three-dimensional Ising model. We already 
know, from section 2, that the interface orthogonal to a lattice axis is rigid 
at low temperatures, and that the corresponding boundary conditions give 
rise to a non translation invariant Gibbs state. It is believed, that at higher 
temperatures, but before reaching the critical temperature T c , the fluctua- 
tions of this interface become unbounded when the volume tends to infinity, 
so that the corresponding Gibbs state in the thermodynamic limit is trans- 
lation invariant. The interface undergoes a roughening phase transition at a 
temperature T = T R <T C . 
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Recalling that the rigid interface may be viewed as a two-dimensional 
polymer system (the system of walls), one might expect that the critical tem- 
perature T^ =2 , of the two-dimensional Ising model, is relevant for the rough- 
ening transition, and that Tr is somewhere near T^ =2 . Indeed, approximate 
methods, used by Weeks et al. (1973), suggest Tr ~ 0.53 T c , a tempera- 
ture slightly higher then T^ =2 . Moreover, van Beijeren (1975) proved, using 
correlation inequalities, that Tr > T^ =2 . 

Since then, however, it appears to be no proof of the fact that Tr < T c , 
i. e., that the roughening transition for the three-dimensional Ising model 
really occurs. 

At present one is able to study rigorously the roughening transition only 
for some simplified models of the microscopic interface. Thus, Frohlich and 
Spencer (1981) have proved this transition for the SOS (solid-on-solid) model. 
Moreover, several restricted SOS models, in which the differences of heights 
at two nearest neighbor sites are restrict to take a few number of values, are 
exactly solvable. They present also a roughening transition (this theory has 
been reviewed by Abraham 1986 and van Beijeren and Nolden 1987). 

From a macroscopic point of view, the roughness of an interface should 
be apparent when considering the shape of the equilibrium crystal associ- 
ated with the system. One knows that a typical equilibrium crystal at low 
temperatures has smooth plane facets linked by rounded edges and corners. 
The area of a particular facet decreases as the temperature is raised and 
the facet finally disappears at a temperature characteristic of its orientation. 
The reader will find information and references on equilibrium crystals in the 
review articles by Abraham (1986), van Beijeren and Nolden (1987), Kotecky 
(1989) and Rottman and Wortis (1984). 

It can be argued that the roughening transition corresponds to the disap- 
pearance of the facet whose orientation is the same as that of the considered 
interface. The exactly solvable SOS models mentioned above, for which the 
function r(n) has been computed, are interesting examples of this behavior 
(this subject has been reviewed by Abraham 1986, Chapter VII, see also 
Kotecky and Miracle-Sole 1986, 1987a, 1987b). This point will be briefly 
discussed in the next two sections, where several results on this subject, con- 
cerning the three-dimensional Ising model, will be reported. See Miracle-Sole 
(1995a) for the proofs and a more detailed discussion. 

4.1 The step free energy 

The step free energy plays an important role in the problem under con- 
sideration. It is defined, using appropriate boundary conditions, as the free 
energy associated with the introduction of a step of height 1 on the interface. 
This quantity can be regarded as an order parameter for the roughening 
transition, analogous, in some sense, to the surface tension in the case of the 
usual phase transitions. Indeed, Bricmont et al. (1982) proved that r step > 
if T < T c d=2 , a result analogous to the above mentioned result by van Bei- 
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jeren (1975). On the other hand r step = if T > T c (see, again, Bricmont et 
al. 1982). 

In order to define the step free energy we consider a parallelepipedic 
box A, defined as in section 2 above, and introduce the (step, m) boundary 
conditions, associated to the unit vectors m = (cos^, sin</>) G R 2 , by 

-(■\ _ f 1 if i > or if z 3 = and %\m\ + z 2 m 2 > C4 1) 

1—1 otherwise 

Then, the step free energy, for a step orthogonal to m (such that m 2 ^ 0), is 
r step (</>) = lim lim lim In ^ tep ' m) (A) (4 2) 

vvv Li^oo^^ooLa^oo /3Li Z(±> n o)(A) V ' 

Clearly, this expression represents the residual free energy due to the consid- 
ered step, per unit length. 

When considering the configurations under the (step, m) boundary con- 
ditions, the step may be viewed as a defect on the rigid interface described 
in section 2. It is, in fact, a long wall going from one side to the other side 
of the box A. A more careful description of it can be obtained as follows. 
At T = 0, the step parallel to the axis (i. e., for m = (0, 1)) is a perfectly 
straight step of height 1. At a low temperature T > 0, some deformations 
appear, connected by straight portions of height 1. The step structure, with 
its probability distribution in the corresponding Gibbs state, can then be 
described as a "gas" of these defects (to be called step-jumps). This system 
can be represented as a polymer system on a one-dimensional lattice. This 
description, somehow similar to the description of the interface of the two- 
dimensional Ising model used by Gallavotti (1972), is valid, in fact, for any 
orientation m of the step. It can be shown that the system of step-jumps, at 
low temperatures, satisfies the properties which are required for applying the 
cluster expansion techniques. In this way, the step structure can be studied. 
Actually, the step-jumps are not independent since the rest of the system 
produces an effective interaction between them. Nevertheless, this interac- 
tion can be treated by means of the low temperature expansion, in terms of 
walls, for the rigid interface, applied to the regions of the interface lying at 
both sides of the step. From this analysis (Miracle-Sole 1995a) one gets the 
following result. 

Theorem 4.1. If the temperature is low enough (i.e., ifT<To, where 
Tq > is a given constant), then the step free energy r step (m), exists in 
the thermodynamic limit, and extends by positive homogeneity to a strictly 
convex function. Moreover, r step (m) can be expressed in terms of an ana- 
lytic function of T, which can be obtained by means of a convergent cluster 
expansion. 

In fact, 

r step (m) = 2(K| + |m 2 |) - (l/(3)((\ mi \ + |m 2 |) lnflmxl + |m 2 |) 

- |mi|ln|mi| - |m 2 |ln|m 2 |) - (1/ /3)<p m (f3) (4.3) 
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where ip m is an analytic function of z = e~ for \z\ < e~ "° . The first two 
terms in this expression, which represent the main contributions for T — > 0, 
come from the ground state of the system under the considered boundary 
conditions. The first term can be recognized as the residual energy of the 
step at zero temperature and, the second term, as — (1//3) times the entropy 
of this ground state. The same two terms occur in the surface tension of 
the two-dimensional Ising model (see Avron et al. 1982 and Dobrushin et al. 
1992). By considering the lowest energy excitations, it can be seen that ip m 
is 0(e -4 ^), and also, that the first term in which this series differs from the 
series associated to the surface tension of the two-dimensional Ising model, 
is 0(e" 12 ^). 

4.2 The shape of a facet 

We have seen, in section 1.5, that the appearance of a facet in the 
equilibrium crystal shape is related, according to the Wulff construction, to 
the existence of a discontinuity in the derivative of the surface tension with 
respect to the orientation. In the case of the three-dimensional Ising model, 
a first result concerning this point, was obtained by Bricmont et al. (1986). 
These authors proved a correlation inequality which establish r step as a lower 
bound to the one-sided derivative dr(9) jdd at 9 = 0+ (here r step = r step (0, 1) 
and t{6) = r(0, sin#, cos0)). Thus r step > implies a kink in r(0) at 9 = 
and, therefore, a facet is expected. 

In fact, r step should be equal to this one-sided derivative. This is reason- 
able, since the increment in surface tension of an interface tilted by an angle 
6, with respect to the surface tension of the rigid interface, can be approx- 
imately identified, for 9 small, with the free energy of a collection of steps. 
The density of the steps being proportional to 6>, the distance between them 
becomes very large when 9 is very small. But, if the interaction between the 
steps can be neglected, the free energy of the whole collection of steps can 
be approximated by the sum of the individual free energies of the steps. 

With the help of the methods described in section 4.1, it is possible to 
study these free energies and to control the approximations involved in the 
problem. The following result can then be derived (see Miracle- Sole 1995a 
for the proof). 

Theorem 4.2. For T < T , we have 

dr(9^)/d9\ e=0+ = T^(<j ) ) (4.4) 

i. e., the step free energy equals the one-sided angular derivative of the surface 
tension. 

It is natural to expect that this equality is true for any T less than 
Tr, and that for T > Tr, both sides in the equality vanish, and thus, the 
disappearance of the facet is involved. However, the condition that the tem- 
perature is low enough is important here. Only when it is fulfilled we have 
the full control on the equilibrium probabilities that is needed in the proofs. 
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The above relation, together with Theorem 1.5, implies that the shape 
of the facet is given by the two-dimensional Wulff construction applied to 
the step free energy r step (m). Namely, 

J r ={xGR 2 :x-m<r step (m)} (4.5) 

where the inequality is assumed for every unit vector m in R 2 . Then, from 
the properties of r step mentioned above, it follows that the facet has a smooth 
boundary without straight segments and, therefore, that the crystal shape 
presents rounded edges and corners. 

Appendix A : Cluster expansions 

Let V be a countable or a finite set, the elements of which will be called 
polymers. Let X C V x V be a reflexive and symmetric relation. We say that 
the elements 71,72 G V are incompatible if (71,72) £ X, and we will also 
write 71 7^ 72. If (71, 72) ^ X we say that the two polymers are compatible. 

For each finite set X C V, we define a(X) to be 1 if any two distinct 
polymers 71,72 G X are compatible, and if some distinct 71,72 G X has 
7i 7^ 72- We assume that a complex valued function 0(7), 7 G V is given, 
and call 0(7) the weight (or activity) of the polymer 7. 

Let a finite set V C V be given. The ensemble of all subsets X C V for 
which a(X) = 1 will be called a "gas" of polymers, while 

z(v) = j2 n c^- 1 ) 

xcv -yex 

will be called the partition function in the volume V of this gas (for X = 
the product in the above expression is interpreted as the number 1). 

We introduce the notion of multi-index as the maps X from V into the 
set of non negative integers, such that 

N(X) = *(7) < 00 (A2) 

■yEV 

and define 

X =II<^) X(7) > X\ =11X^)1 (A3) 

We denote by suppX the set of polymers for which -^(7) 7^ 0. A multi-index 
can be viewed as a finite family of elements of V, X(^j) being the multiplicity 
of the element 7 in the family X. 

To any multi-index one associates a graph G(X) with N(X) points: 
X(7i) distinct vertices associated to 71, etc.. We put a line between 71,72, 
if, and only if, 71 and 72 are not compatible (in particular if 71 = 72 we put 
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a line between them). A subgraph C of G(X) is called a full subgraph if the 
sets of vertices of C and of G(X) coincide. 

Take a T (X) = if G(X) is disconnected, and 

a T (X) = (X!)- 1 £ (-!)*<*> (A5) 

CcG(X) 
C connected 

otherwise. Here the sum runs over all connected full subgraphs of G(X) and 
£(X) is the number of lines of C (if G(x) has a single point we interpret the 
sum as having one graph with £(C) = 1). 

Theorem A.l. We have 

lnZ(V)= CL T {X)(j) X (AA) 

X: suppXcV 

The connected X will be called clusters (of polymers). The expansions 
in terms of these objects are the cluster expansions. Their convergence prop- 
erties will be established in the Theorem below. If the weights 4>{l) are small 
enough for "large" polymers 7, this Theorem enables one to make good eval- 
uations of the free energy as well as of the decay of correlations. Let us 
remark also that whenever is analytic in some parameter the estimates 
stated below imply the analyticity of the free energy and the correlation 
functions. 

Theorem A. 2. Assume that there exists a positive real valued function 
^(7), 7 G V , such that 

|0( 7 )| ^(7) _1 exp ( M7)) < r < 1 (A6) 

for each 7 G V . Then the following estimate 

|a T (X)0 x |<M7)r(l-r)- 1 {A.l) 

X: 76suppX 

holds true for every 7 G V . 

See Gallavotti et al. (1973), chapter 4, for the proof of the results stated 
in Theorems A.l and A. 2, in the particular case of the Ising model at low 
temperatures. In this case the term contour is used instead of polymer. 
To stress that the formulation, as well as the proof of these results, do not 
depend on details of the "geometry" of polymers, we have presented them in 
an abstract setting. 
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2.3 Appendix B : Proof of Theorem 2.1 

In this appendix we prove Theorem 2.1 for the three-dimensional Ising 
model. We write 

F(L lt L 2 , L 3 ) = F(A) = In (^±> n ) (A)/Z+(A)) (fl.l) 

First, we establish the monotonicity property 

F(L U L 2 ,L 3 ) <F(Lx,L 2 ,4) if L 3 > 4 (S.2) 

for Z7 d large enough. This shows that the limit L 3 — > oo exists and 

lim F{L U L 2 , L 3 ) =inf F(L U L 2 ,L 3 ) = F(L U L 2 ) (B.S) 

We shall then prove that F(Li, L 2 ) is bounded 

0<F(Lx,L 2 ) <K X L X L 2 (BA) 

and that it satisfies the following subadditivity property (up to boundary 
terms) 

F(L[ + L'l + 1, L 2 ) < F(L;, L 2 ) + F(L'{, L 2 ) + K 2 {L X + L 2 ) (B.5) 

The subadditivity (B.5) together with the bound (B.4) imply, following stan- 
dard arguments in the theory of thermodynamic limits (see for instance Ru- 
elle 1969) the existence of 

lim J—F{L U L 2 ) (B.6) 

Li ,iv2— >-cxd Li\1j 2 

Moreover this limit equals the infimum over Li,L 2 . This ends the proof of 
the first part of Theorem 2.1, provided properties (B.2), (B.4) and (B.5) are 
satisfied. 

The proof of these properties uses Griffiths correlation inequalities (see 
for instance Ruelle 1969), and is based in the following Lemma. This same 
Lemma will be used to prove the second part of Theorem 2.1. 

Lemma. Given two boxes A, A' such that A' C A we let S\ be the set 
of sites of A \ A' which are above or on the the plane p n and S 2 the sites of 
A \ A' which are below this plane. Then, we have 

< F(A) < F(A') + 2/3N(Si, S 2 ) (B.7) 

where N(Si, S 2 ) is the number of bonds with i e S\ and j E S 2 . 
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The monotonicity property (B.2) is just a particular case of the above 
inequality (B.7). If A and A' are the two boxes considered in (B.2) one has 



The upper bound (B.4) follows also from inequality (B.7) by taking A' = 
and noticing that in this case N(Si, S2) < 3LiL 2 , which gives K = 3/3 J. 

In order to prove the subadditivity property (B.5) we consider three 
parallelepipeds Ai, A2 and A of sides (L' 1 ,L2,L 3 ),(L'{,L2,L 3 ) and (L[ + 
L'l + 1, L2, L3) placed in such a way that Ai U A2 C A, the distance from Ai 
to A2 is 1, and the plane p n which passes through the center of A passes also 
through the center of Ai and A2, and defines the (±) boundary condition 
for the three boxes. Now we apply inequality (B.7) with A' = Ai U A 2 . We 
notice that A\A' reduces to the vertical plane separating Ai and A2 when the 
vertical dimensions L' 3 , L 3 and L 3 , tend to infinity. Thus Af(Si, S2) < 2Li, 
and F(A'), the first term in the right hand side of (B.7), nearly coincides 
with the sum F(L / 1 ,L 2 ) + F(L", L 2 ) because the partition function Z tT (A') 
factorizes. 

Some error is made because from the construction the centers of Ai and 
A2 do not necessarily coincide with a site of the lattice, as it is assumed in 
the definition of r(n). This error however is bounded by the third term in 
the left hand side of (B.5) as it may be seen by using appropriate planes 
parallel to p n introduced to compensate the displacement of the center (by 
a distance less than one) and applying again the above argument. 

We next prove that r(n) satisfies the pyramidal inequality. 

Let Aq, . . . , A 3 be the vertices of the pyramid. Introduce their projec- 
tions A'q,..., A' 3 on the horizontal plane, and assume that A' falls in the 
interior of the triangle A[, A' 2 , A' 3 . Let C! be the projection of the lattice C. 
We denote by Qq the set of sites of CJ inside the triangle A[, A 2 , A 3 and by 
Qi the set of sites of C inside the triangle A' ,A' 2 ,A 3 whose distance from 
the sides of this triangle is larger than 1. Similarly we define the sets Q2 and 
Q 3 with respect to the triangles A^A'^A^ and A^A^A 2 . We introduce 
the surface which coincides with the plane Ai,A 2 ,A 3 outside the triangle 
Ai, A 2 , A 3 and with the other three faces of the pyramid inside it. We apply 
now a modified version of the Lemma in which this surface replaces the plane 
p n . Assume that Q is the projection of the box A and define A' as the set 
of sites in A whose projections belong to Qi U Q 2 U Q3. Then the modified 
Lemma implies 



F(Q ) < F{Qi) + F(Q 2 ) + F(Q 3 ) + \Q \ {Q 1 UQ 2 U Q 3 )\ (K/a) (B.8) 



From this relation, the pyramidal inequality follows by passing to the limit 
when the three triangles tend to infinity. 

Proof of the Lemma. To prove inequalities (B.7) we notice that F(A') 
may be obtained from F(A) by adding the external fields 



N(S 1 ,S 2 ) = 0. 




ieSius 2 
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4^(r(j) + /i^(T(i) to H(a A \±) (B.9) 

ieSi ies 2 

and letting h tend to infinity. In fact after this limit the partition functions 
Z"(A) in (B.l) are replaced by Z^(A') multiplied by the term 

exp ( Yl ^W^O')) 

(i,j)cS 1 US 2 

This means that F(A) becomes equal to the right hand side of equation (B.7). 
But, Griffiths inequalities tell us that 

(a(z))+-|(a(z)) ± |>0 (B.10) 

where ( ) a denotes the expectation values corresponding to the Gibbs mea- 
sure (Z ") -1 exp(— /3H(a \ a), and show that the derivative with respect to h 
of the modified F(A) is positive, and hence that it is an increasing function 
of h. We get therefore the second inequality in (B.7). The same argument 
implies that Z + (A) is greater than Z ± (A) and gives the first inequality in 
(B.7). 
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